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Abstract. The formation of mounded surfaces in epitaxial growth is attributed 
to the presence of barriers against interlayer diffusion in the terrace edges, known 
as Ehrlich-Schwoebel (ES) barriers. We investigate a model for epitaxial growth 
using a ES barrier explicitly dependent on the step height. Our model has an 
intrinsic topological step barrier even in the absence of an explicit ES barrier. 
We show that mounded morphologies can be obtained even for a small barrier 
while a self-affine growth, consistent with the Villain-Lai-Das Sarma equation, 
is observed in absence of an explicit step barrier. The mounded surfaces are 
described by a super-roughness dynamical scaling characterized by locally smooth 
(faceted) surfaces and a global roughness exponent a > 1. The thin film limit 
is featured by surfaces with self- assembled three-dimensional structures having 
an aspect ratio (height /width) that may increase or decrease with temperature 
depending on the strength of step barrier. 
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1. Introduction 

The most of the experiments of film production is performed in conditions far from 
equihbrium, in which kinetic phenomena usuaUy overcome any thermodynamic effects 
in the film properties. Morphology, one of the most important features in the 
film production, is ruled by distinct diffusion mechanisms [1, 2], deposition rate [3] 
and temperature [4, 5] among other factors. In particular, the emergence of self- 
assembled and/or three-dimensional structures in form of mounds has been observed 
during the growth of a wide diversity of films, ranging from metals [6, 7, 8] to 
inorganic [9, 5] and organic semiconductor materials [10, 11]. This mound instability 
may be an undesirable feature in the production of atomically fiat interfaces but 
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becomes a valuable tool to produce structured surfaces, both systems with important 
technological applications. 

The three-dimensional growth is attributed to kinetic barriers that promote the 
unbalance between up and downhill currents in stepped surfaces. A number of different 
mechanisms has been proposed to explain these barriers, but there are only few cases in 
which the observation of an unstable growth can clearly be associated to a particular 
barrier [5]. In two seminal papers, Ehrlich et al. [12] and Schwoebel et al. [13] 
independently reported that adatoms lying in the edge of terraces diffuse backwardly 
in the same layer more frequently than move to a lower step. This experiment 
led to a widely accepted hypothesis of an additional activation energy to downhill 
diffusion known as Ehrlich- Schwoebel (ES) barrier. Besides the ES barrier, alternative 
mechanisms have been proposed to explain mound instability. Examples include the 
short-range attraction of adatoms towards ascending steps [14] and fast edge diffusion 
[15, 16, 17]. 

Kinetic Monte Carlo (KMC) is a standard simulation method used to investigate 
morphology and dynamics during the epitaxial film growth [2] . In a KMC simulation, 
the events (deposition, diffusion, evaporation, etc.) are implemented, one at a time, 
with the specific rates. In particular, diffusion is a thermally activated process 
occurring at a rate given by an Arrhenius law D = z/q exp(— E'l^/Zc^T), where z/q 
is an attempt frequency, T the growth temperature and Ed the diffusion activation 
energy [18]. An ES barrier is usually included in KMC simulations by raising Ed 
when the particles are lying in descending steps. This additional energy effectively 
inhibits the downhill diffusion and, consequently, increases the chance of nucleation 
of new layers in top of the terraces. The result is an instability responsible by the 
emergence of three-dimensional structures [19]. 

Molecular dynamic simulations show that an efficient upward diffusion can also 
be present in epitaxial growth [20]. Additionally, molecular dynamic simulations also 
show that the ES barrier can change substantially from monolayer to multilayer steps 
[21]. With these motivations, we investigate a model with a step-height dependent 
ES barrier using KMC simulations. We are interested in the effect of these barriers in 
both down and uphill diffusion in multilayer steps. We have organized the paper as 
follows. In section 2, we develop the model and present the simulation procedures. In 
section 3, the simulation results are presented and discussed. Finally, our concluding 
remarks are done in section 4. 

2. Model 

The KMC simulations were performed as follows. Particles are deposited on a 
substrate kept at a constant temperature T. The substrate, represented by a triangular 
lattice with periodic boundary conditions, is initially fiat. Deposition occurs normally 
to the substrate at a constant rate F and obeys the solid-on-solid restriction to prevent 
overhangs. The height profile hj is therefore the number of particles adsorbed in the 
site j. The deposition involves two steps: firstly a site is picked up at random and 
secondly a new particle is deposited in the most energetically favourable (largest bond 
number) site among the chosen one and its nearest neighbours (NN). If there are 
multiple options, one is randomly chosen. This deposition rule is equivalent to the 
classical Wolf- Villain model [22]. 

Intralayer diffusion from a site j to j' ^ Ah{j^j') = hj' -h 1 — /ij = 0, has an 
activation energy ED{j^f) = Eq -\- njEN, where Eq represents the interaction with 
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Figure 1. The barrier model in multilayer steps, (a) A particle in an ascending 
step detaches from its initial position (empty square) and starts an unbiased 
random walk along the normal direction (red), (b) The particle can be reflected 
in the top due the presence of a ES barrier, (c) A particle in a descending step 
(empty square) overcomes the ES barrier with probability p = exp(— ^^b/fc^T) 
and then diffuses exactly as in (b). (d) Random walk problem including the 
absorbing/reflecting boundaries indicated by X. 



the substrate and the contribution of each one of the nj lateral bonds. The 
diffusion rate is therefore given by the Arrhenius law: 



(1) 



Apart of the step barrier, this diffusion rule is similar to that used in the model 
proposed by Smilauer and Vvedenski [18]. Interlayer diffusion is allowed for particles 
lying in both ascending or descending steps. For monosteps, Ah{j^j') = ±1, the 
barrier is implemented as usually: an additional barrier £^5 is present in a descending 
but not in an ascending step. For steps higher than a monolayer, the interlayer 
diffusion is implemented in two parts. If the diffusing particle lies in an ascending step, 
it firstly detaches from its initial layer at a rate given by Eq. (1) and then executes 
a non-directed one-dimensional random walk, perpendicularly to the substrate, as 
schematically illustrated in Fig. 1(a). A step barrier is also present in the top of the 
step [Fig. 1(b)]. The rule for a particle in a descending step is very similar except 
that it has to overcome the ES barrier before attaching to the lateral [Fig. 1(c)]. 
Notice that the barrier in the top makes the model symmetric in relation to down and 
upward diffusion in multilayer steps, an improvement of the usual rule without height 
restrictions. 

The determination of the probability of the particle moving to a neighbour site or 
being reflected to its initial position involves the solution of a one-dimensional random 
walk with absorbing/reflecting boundaries. When the particle visits the boundaries, 
it becomes trapped with a given probability or is reflected with complementary 
probability. This problem is schematically illustrated in Fig. 1(d) for a step of height 
i = 4. In our model, each boundary represents either the initial or neighbour sites. 
The reflection is due to the ES barrier and, therefore, we set adsorption probabilities 
1 and p = exp(— ^^^//c^T) for the the bottom {i = 0) and the top (i = + 1) of the 
step, respectively. Let q{£ -\- l\i) be the probability of the particle be trapped in the 
top if it starts at a distance i from the bottom. This random walk problem has an 
exact solution explicitly developed in Ref.[23]: 

.(^+1W = ^. (2) 

For an ascending step [Fig. 1(a)], it is easy to see that the probability of moving 
to the other layer is q{i -\- 1|1), while for a descending step this probability becomes 
q{£ X p. In latter case, the factor p takes into account the step barrier when the 
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Figure 2. (Colour on-line) Surface morphologies for growth times varying from 
5 to 400 ML. The growth temperature is T = 573 K and a mild step barrier 
Eb = 0.02 eV is used. 

particle crosses the kink [Fig. 1(c)]. The resulting probability of diffusion through the 
step, for both cases, reads as: 

= 1+,(|L|-1) ' (3) 
where Ah = hj' -\- 1 — hj is the height of the step to be crossed. 

3. Results 

The KMC simulations were performed with fixed parameters = 1 eV, = 
0.11 eV, F = 1 monolayers (ML) per second, = 10^^ s~^. Temperatures 
T = 473, 523 and 573 K were simulated for deposition times of 10^ ML. Except when 
explicitly mentioned, we report results for 256 x 256 sites. The initial surface evolution 
for a growth temperature of T = 573 K and a weak ES barrier of £^5 = 0.02 eV can 
be followed in Fig. 2. One can see the surface coarsening and the emergence of three- 
dimensional structures with a characteristic size in the form of mounds. For very 
long times, however, the instability is highly enhanced turning the surfaces in deeply 
grooved morphologies with large terraces. 




(a) (b) 

Figure 3. (Colour on-line) (a) Surface morphology for distinct growth 
temperatures after the deposition of 100 ML. The bottom and top surfaces 
correspond to step barriers Ei^ = 0.02 eV and Eij = 0.07 eV, respectively, 
(b) Grain morphology for E^ = 0.07 eV and growth temperatures T = 523 K 
(top) and T = 573 K (bottom). Both grains are in the same scale. 
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Figure 4. (Colour on-line) Characteristic length against time for (a) T = 473 K 
and (b) T = 523 K using distinct step barriers. Inset shows the correlation 
functions for different deposition times and = 0.07 eV. The curves correspond 
to averages over 20 independent samples. The error bars smaller than symbols 
were omitted. Dashed lines are power laws with the indicated exponent as guides 
to the eyes. 



A rich variety of spatio-temporal properties is observed in the thin film phase. 
The present model has an intrinsic step barrier independent of the parameter E^j when 
\Ah\ > 2. If we set ^5 = 0, the system has a transition from a rough (fractal) surface at 
low temperatures to an almost layer- by-layer growth at high temperatures [inset of Fig. 
5(b)]. The layer- by- layer regime is transient and becomes self-affine for asymptotic 
times, as we will see later. However, a weak barrier of Eij = 0.02 eV is sufficient to 
induce a temperature-driven instability and to form mounds. In the top row of Fig. 
3(a), we can see the temperature effect for this weak barrier. The mound aspect ratio 
(height/width) decreases with temperature, a behaviour also observed in homoepitaxy 
of some semiconductor materials [24]. If a step barrier ^5 = 0.07 eV is used instead, 
the surface morphology changes considerably. The uphill diffusion becomes strong at 
higher temperatures generating self-assembled three-dimensional structures, as shown 
in the bottom of Fig. 3(a). Details are illustrated in Fig. 3(b), where isolated grains 
for temperatures T = 523 and 573 K are shown. Higher temperatures make the grains 
higher resulting a large increase of the aspect ratio. 

The thick film limit also exhibits very interesting dynamics. A basic quantity to 
characterize the surface morphology and dynamics in this regime is the height-height 
correlation function defined as [15]: 

r(r) = {/i(x)Mx + r))x, (4) 

where h{x.) is the surface height in the mean height referential and (• • •)x represents the 
average over the surface. The first zero of r(r), denoted by ^, is a characteristic length 
of the surface. Typical correlation functions for distinct growth times are shown in the 
insertion to Fig. 4(a). The oscillating behaviour observed in the correlation functions 
is the hallmark of mounded surfaces [1, 15]. The length ^ is expected to increase as 
a power law of time, ^ ~ t*^, where ( is the so-called coarsening exponent [15]. Other 
basic quantity to characterize the surface dynamics is the interface rms roughness 
given hy w = ^/T(fi). For a kinetic roughening, this quantity grows as a power law 
of time, w ^ t^ , where f3 is the growth exponent. In dynamical scaling theory [1], 
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Figure 5. Interface rms roughness evolution for (a) T = 473 K and (b) 
T = 523 K. The inset shows the initial layer-by-layer growth for T = 523 K. 
Dashed lines are scaling laws obtained from data regressions. 



the coarsening and growth exponents are related with the roughness exponent by 
the scahng relation a = /3/(. These exponents are very useful because they may 
determine the universality class of the system and, from the framework of the generic 
dynamical scaling theory, they may connect scaling and morphological properties of 
the surface [25]. 

For all investigated temperatures the scaling exponents change if an explicit step 
barrier is present or not. For the lower temperature of T = 473 K and in the absence 
of an explicit step barrier {E^ = 0), we found a power law regime with coarsening and 
growth exponents ( = 0.32 and f3 = 0.19, respectively, as shown in Figs. 4(a) and 5(a). 
A surface evolution for these parameters is shown in MOVIE 1 of the supplementary 
material. These exponent values are in good agreement with the exponents ( = 3/10 
and P = 1/b obtained for the classical Villain-Lai-Das Sarma (VLDS) growth equation 
[19, 26] 

^ = -vV^h + AV(V/i)2 + r?, (5) 

where 77 is a Gaussian noise. Similar exponents are obtained for higher temperatures 
and El) = 0. Interestingly, in the limit of zero temperature, our model corresponds 
to the Wolf- Villain model [19] that has a crossover [27, 28] to the universality class 
of the Edwards- Wilkinson (EW) equation {dth = vV'^h + rf) with exponents /3 = 
(logarithmic growth) and C = 0-5in(i = 2 + l [29] . Contrarily, our simulations 
at higher temperatures point out the opposite behaviour: a small growth exponent 
(consistent with a logarithm growth) and a coarsening exponent close to 1/2 are 
observed for early growth times whereas /3 1/5 and C ~ 3/10 are seen only at long 
times. The results suggest that the diffusion, not the deposition rules, is governing 
the dynamics. 

In the presence of an explicit ES barrier, the scaling exponents for T = 473 K and 
T = 523 K converge (more slowly for Ei^ = 0.02 eV than for 0.07 eV) to (3 = 0.33 and 
C = 0.22, as shown in Figs. 4 and 5. A growth exponent /3 = 0.33 were also obtained 
for the higher investigated temperature of T = 573 K, but the measured coarsening 
exponent, ( ~ 0.18, seems not be the asymptotic value. Anyway, the corresponding 
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Figure 6. (Colour on-line) Surface morphologies after deposition of 10^ ML using 
an step barrier = 0.07 eV at temperatures T = 473 K (left) and T = 523 K 
(right). 

roughness exponents, given by a = are larger than 1 implying that the explicit 
barrier causes the so-called super-roughness dynamical regime [25]. Super-roughness 
is featured by a locally smooth surface with rms roughness in windows of size £ scaling 
as w{£) ^ £ for £ « ^. Left panel of Fig. 6 shows surfaces for a long deposition 
time obtained for = 0.07 eV. Similar structures are obtained for a mild barrier 
E}) = 0.02 eV, but the faceted morphology is less evident. For T — 473 K, the 
system evolves from a initially rough [Fig. 3(a)] to a faceted morphology exhibiting 
pyramidal mounds as can be seen in MOVIE 2 of the supplementary material. The 
surface dynamics becomes richer for higher temperatures. Initially, self-assemble and 
fast growing structures are formed due to the instability caused by the step barrier. 
Thereafter, these structures start to coalesce and to form large scale mounds also 
exhibiting pyramidal morphology. The surface evolution can be followed in MOVIE 3 
of the supplementary material. 

4. Concluding remarks 

The pattern formation in semiconductor films produced with epitaxial techniques is 
a subject of basic and technological interest. The formation of three-dimensional 
structures is related with the so called Ehrlich-Schwoebel (ES) barrier. This barrier 
is usually implemented in kinetic Monte Carlo (KMC) simulations by reducing the 
diffusion rate of particles in descending steps. In the present work, we investigate 
the homoepitaxial growth of films performing KMC simulations where the ES barrier 
depends explicitly on the step height. In addition, our model includes a barrier to 
inter layer diffusion in ascending steps when their heights are larger than a monolayer. 

The model yields surfaces with mounded morphology even for a mild step barrier 
of E\) = 0.02 eV, a value much smaller than the typical nearest neighbour interactions 
of E^ = 0.11 eV used in the present KMC simulations. However, a layer-layer growth, 
lasting during the deposition of several layers, is obtained if the explicit ES barrier 
is null. For an intermediary barrier E\) = 0.07 eV, self-assembled three-dimensional 
structures are obtained. The aspect ratio (height /width) of these structures decreases 
with temperature for a mild barrier but increases for the intermediary one. 

We have also investigated the thick film limit using the growth and coarsening 
exponents, defined as w ^ and ^ t^, respectively. In the absence of an 
explicit barrier, the interface undergoes kinetic roughening with exponents /3 ~ 0.19 
and C ~ 0.30 what is consistent with the non-linear Villain-Lai-Das Sarma (VLDS) 
equation [19, 26] for molecular beam epitaxy (MBE). In the presence of a step barrier. 
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the exponent changes to 13 ^ 0.33 and C ~ 0.22 what turns out a roughness exponent 
a = /3/( > 1 and states the so-cahed super-roughness scahng regime [25] for the 
surface evolution. This regime is featured by large pyramidal structures that are 
observed in MBE systems like Ge(lOO) [24] and Cu(lOO) [30]. Even though a large 
number of models has been proposed to reproduce specific systems with pyramidal 
morphology [2], our model is quite simple, robust and does not require any kind of 
parameter tuning. 

The present model addresses two features that, to our knowledge, were neglected 
in the broad literature of KMC modelling of epitaxial growth: the influence of the step 
height in the ES barrier and the explicit barrier when the particle is executing an uphill 
diffusion. This last one is particularly interesting since the rules for the movement 
through a kink becomes symmetric in relation to up or downhill movements [Fig. 
1(b)]. Finally, the specific form of the interlayer diffusion probability [Eq. (3)] is not 
essential to the surface dynamics. We have also simulated, but not shown in this 
paper, a model where Eq. (3) is replaced by Pcross = P for |^^| ^ 2 and observed 
that all morphological properties and scaling exponents are preserved. 
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Supplementary material of "Modelling of epitaxial 
film growth with a Ehrlich-Schwoebel barrier 
dependent on the step height" 
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Abstract. We described the animations of the surface evolution generated with 
kinetic Monte Carlo simulations. 



MOVIE 2 T473Eb00L256_sinall . avi - Surface evolution for a substrate temperature 
of T = 473 K in the absence of an explicit step barrier. 

MOVIE 2 T473Eb07L256_sinall . avi - Surface evolution for a substrate temperature 
of T = 473 K and a step barrier of = 0.07 eV. 

MOVIE 3 T523Eb07L256_small . avi - Surface evolution for a substrate temperature 
of T = 523 K and a step barrier of Ei, = 0.07 eV. Scale changes during the movie 
for sake of visibility. 



